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Abstract 



Q , In this paper, we study in details the critical behavior of the 0{n) quantum tp^ 

^ I model with long-range interaction decaying with the distances r by a power 

law as r~ in the large n-limit. The zero-temperature critical behavior 

is discussed. Its alteration by the finite temperature and/or finite sizes in 

^ • the space is studied. The scaling behaviours are studied in different regimes 

depending upon whether the finite temperature or the finite sizes of the system 

is leading. A number of results for the correlation length, critical amplitudes 

and the finite-size shift, for different dimensionalities between the lower (i< = 

'^ ' a/2 and the upper (i> = 3a/2 critical dimensions, are calculated. 
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I. INTRODUCTION 

The vast majority of existing analytical and numerical work on finite-size effects has been 
in regimes where the quantum effects are entirely ruled out [|l|]. In this case both the statics 
and the dynamics can be described by classical statistical models. At the critical point the 
bulk thermodynamic functions of these models are singular. For these systems the finite-size 
scaling theory asserts that the singularities holding at the thermodynamic limit are altered 
depending upon the nature of the geometry to which the system is confined and the imposed 
boundary conditions [0,§ . Indeed various type of geometries, can be considered, depending 
on the number of the finite sizes in the model. 

The 0{n) vector models are extensively used to explore the finite-size scaling theory, 
using various methods and techniques. For finite n the most used method is that of renor- 
malization group |^. The most investigated case is the one corresponding to the limit 
n = oo (this limit includes also the mean spherical model) |^. In this limit, these models 
are exactly soluble for arbitrary dimensions and in a general geometry. The most part of 
these investigation are limited to systems in which the forces are short-range. To test the 
finite-size scaling idea, when the interaction is of a long range (varying with a power law), 
the only used model is the mean spherical model f^. 

In recent years there has been an increased interest in the theory of zero temperature 
quantum phase transitions P,0]. Distinctively from temperature driven critical phenomena, 
these phase transitions occur at zero temperature as a function of some non-thermal control 
parameter (or a competition between different parameters describing the basic interaction of 
the system) , and the relevant fluctuations are of quantum rather than thermal nature. In this 
type of critical phenomena the time plays a crucial and fundamental role. The coupling of 
statistics and dynamics that is inherent to quantum statistical problems introduces effective 
dimensionalities in the hyperscaling laws, i.e. the space dimensionality d is replaced by ci-l-^ 
[z is the dynamic critical exponent). In this case the inverse temperature acts as a flnite 
size in the "imaginary-time" direction for the quantum system at its critical point. This 
allows the investigation of scaling laws for quantum systems near the quantum critical point 
in terms of the theory of flnite-size scaling [Q or using the conformal fleld theory techniques 
by mapping the bulk system in a flnite one [^. The 0{n) symmetric vector models are also 
used in the exploration of of the quantum critical phenomena and for the investigation of the 
scaling properties of such phase transitions. The quantization of classical 0{n) is performed 
with the help of the Trotter formula which maps the quantum model on a classical one with 
z additional effective dimensions. 

In systems showing quantum critical behaviour the temperature plays two different roles. 
For temperatures low enough, quantum effects are essential. In this case the temperature 
affects the geometry to which the system is conflned adding a "new" sizes to the eucludean 
space-time coordinate system. By raising the temperature, the system is driven away from 
quantum criticality. At high temperatures, however, the size in the "imaginary-time" direc- 
tion becomes irrelevant in comparison with all length scales in the system. In this case we 
have a classical system in d dimensions and the temperature is just a coupling constant in 
the classical critical behaviour. 

In this paper we present a detailed investigation of the scaling properties of the quantum 
0{n) vector ip"^ model with long-range interaction. Our study will include the quantum as 



well as the finite-size effects and their influence on the critical behaviour and the critical am- 
plitudes. We will also check the influence of the interaction range on the critical behaviour. 
These interactions enter the exact expressions for the free energy only through their Fourier 
transform which leading asymptotic is U{q) ~ q'^* , where a* = min{a,2) |jlO|. As it was 



shown for bulk systems by renormalization group arguments a > 2 corresponds to the case 



of finite (short) range interactions, i.e. the universality class then does not depend on cr [jlO 
Values satisfying < cr < 2 correspond to long-range interactions and the critical behaviour 
depends on a. On the above reasoning one usually considers the case a > 2 as uninteresting 
for critical effects, even for the finite-size treatments [O . So, here we will consider only the 
case < cr < 2. 

The paper is devised as follows: in Section ||, we present some predictions which extend 
the finite-size scaling to quantum systems. We also discuss the interplay of quantum and 
finite-size effects on the quantum critical behaviour. We comment on the anisotropy caused 
by the presence of both effects. In Section |T|, we review, briefly, the (p^-model with long- 
range interaction and present the saddle point equation. We investigate the low temperature 
behaviour of the bulk model in Section |V|. In Section [V|, we investigate in details the finite- 
size behaviour at zero temperature. Section ^ is devoted to some comments about the low 
temperature and finite-size effects of the system. In Section |Vli] we discuss our result briefly. 
In the remainder of the paper we present some details of the calculations. 

II. FINITE-SIZE SCALING AND QUANTUM CRITICAL BEHAVIOUR 

Divergent length scales play a crucial role in continuous phase transitions. Unlike classical 
models, where the scaling can be done equally for all "spatial" dimensions, quantum models 
are anisotropic in general, and therefore the "space" and "imaginary-time" directions will not 
scale in the same fashion. According to general hypothesis of finite-size scaling theory |T^ 
extended here for quantum (anisotropic) system, a physical quantity A (r, h, L, T) (with r 
measuring the distance from the critical point, h is an ordering field coupled to the order 
parameter, T is the temperature and L is the size of the system), which may be singular at 
the critical point in the thermodynamic (bulk) limit at the quantum critical point (r = 0), 
will scale like 

A{r, h, T, L) = FAs (r6^/^ /^6^/^ T6^ ftL"^) . (2.1) 

In the scaling form ( ^.l] ), p corresponds to the engineering dimension d+ z oi the system in 
the case when the scaling function refers to the singular part of the free energy, and it is the 
divided by u critical exponent measuring the divergence of the bulk thermodynamic function 
A at the critical point for the other physical quantities of interest (for the correlation length 
p = 1, for the susceptibility p = 7/z/, etc). Depending upon the choice of the renormalization 
group rescaling factor b we obtain different scaling functions As, which are related among 
each other by some appropriate change of the scaling variables. 

Before starting to discuss the general form of the scaling universal function As, we will 
discuss the two limiting cases which were subjects of several investigations during the past 
two decades. 



The first case corresponds to zero temperature and is called hereafter the quantum critical 
behaviour. Here equation ( p.l|) reduces to 



A{r, h, 0, L) = IfA^ irh^^", hb^/", bL^^) . (2.2) 



Following Ref. [13], we choose the rescaling factor b to be proportional to the linear size of 



the system. Then, we obtain 

Air, h, 0, L) = L^A^ (rL^'\ hL^'^) . (2.3) 

Here the situation resembles that of systems exhibiting classical (thermal) phase transition. 
It is possible to get another result for the scaling function in the right hand side of 
equation ( |2.3| ), if one considers the variable fL^/'^ instead of the first variable of the scaling 
function A^. The parameter f is introduced in such a way to account for the shift of the 
critical quantum parameter to the value corresponding to the rounding of the thermody- 
namic functions, when the number of the infinite dimensions d' is less than its lower critical 
dimension d<. For the opposite case we have just a shift of the critical control parameter. 
In this case we find that the critical exponents are those of a d' dimensional system (see for 
example |@|). By definition the finite-size shift is given by 

f = r + e(L), lime(L) = 0. (2.4) 

L— *oo 

In general we have 

e{L) ~ L-l/^ (2.5) 

and so this quantity shrinks to zero in the thermodynamic limit. When the arguments 
of the scaling functions get replaced by zero, we will obtain universal critical amplitudes, 
characterizing the whole class of universality. Here we have to emphasize that the values of 
the universal critical amplitudes are different depending upon the point in which they are 
calculated i.e. r or f. 

The second case we consider is the one corresponding to the bulk system (L = oo) at 
finite-temperature. In this case the scaling form ( |2.1| ) transforms into (see Ref. ^ and 
references there) 

A{r, h, T, oo) = T~P/'Al (rT-^/"", hT-^l'^) , (2.6) 



where we used the relation b = T^ between the temperature and the rescaling parameter b. 
The same predictions for the shift and the critical amplitudes remain valid here. We find it 
convenient to use Lr ~ T~^ as a linear "temporal" size instead of the inverse temperature. 
The general case corresponding to both quantum and finite-size effects can be studied by 
considering a system with larger dimension p = d + z confined to a general geometry of the 
form L'P~^ X L^ WM . The fact that the inverse temperature can be used as an additional size in 



the imaginary-time direction creates some anisotropy in the system. This property will lead 
to the establishment of some change in the scaling properties of the finite quantum system. 
In the general case we can consider the quantum to classical and the finite-size to the bulk 



system different crossover phenomena. It is easy to convince oneself about this statement by 
thinking on the dynamic critical exponent which is different for different quantum systems. 
The situation of the combined investigations of finite-size scaling and quantum to classical 
crossover is similar to the one formulated in the framework of the phenomenological study of 
finite-size scaling in anisotropic systems [jl5|. In this case one can investigate the interplay 
of quantum and finite-size scaling. For example in the case when the quantum effects are 
leading compared to the finite-size (this case will be called here "low-temperature" case) i.e. 
Lt- <C L, from equation ( |2.1| ) to first order in L^/L, we expect to obtain 



+Ll+^L-^Al'' {rLy\ hL^/") (2.7) 

instead of equation ( |2.6| ). This shows how the finite-size effects give rise to some corrections 
to the quantum scaling. Following the same reasoning in the case of "very low-temperature" 
i.e. Lr ^ L, when the finite-size contributions to the ground state energy and its derivatives 
are dominant compared to those coming from the quantum effects, we find 

^(r, /., T, L) = LMf (rLy^ /iLf Z'^) 

+L;'LP+'A^^ {rL\'\ hL^'^) , (2.8) 

showing a correction to the zero-temperature finite-size scaling due to the temperature. 
These ideas has been tested in the framework of the spherical quantum rotor model in 
Ref . [0,11^ . We will compare the scaling forms (p.7|) and (|2.8| ) with the available analytical 



results obtained below. 



III. THE FREE ENERGY AND THE GAP EQUATION 

The quantum yj^-model we are going to investigate its quantum critical and finite-size 
scaling properties is (for a review of the applicability of this model in exploring quantum 
critical phenomena see Ref. R]) 



'H{^} = \ fj^ dr J dx [{drvf + (V-/V) ' + W + ^V' 



(3.1) 



where </? is a short hand notation for the space-time dependent n-component field y9(x,r), 
uq and tq are model constants and T is the temperature. In (|3.1|) we assumed h = ks = 1 
and the size scale is measured in units in which the velocity of excitations c = 1. Here we 
will consider periodic boundary conditions. This means 

(p{x, t) = \ —Y^ ip{k, uji) exp {tk- X - zujit) , (3.2) 

where oji = 2ttTI (with / = 0, ±1, ±2, ■ • •) are the Matsubara frequencies for bosonic systems, 
A; is a discrete vector with components ki = 27rnj/L, n, = 0, ±1, ±2, ■ ■ ■ and a cutoff A, and 
V = L'^ is the volume of the system. We note that the second term in the model transforms 



into \k\"(p^{k^ uo) in the momentum representation, where the parameter < a < 2 accounts 
for short-range and long-range interaction as well. 

The partition function of the Hamiltonian (p.l|) reads 



(3.3) 



Z= f V(p exp {-n {ifi}) . 



Note that in the low-temperature limit T -C A the integral over r in 7^{0} (see (pTl|)) can be 
extended over the whole temperature axis to give an effective (p"^ model in d + z dimensions 
with a quantum control parameter tq. In the high temperature limit T 3> A the upper limit 
in the integral over r is small. This offers us the possibility to write the Hamiltonian as a 
classical ip"^ model in d dimensions. Now using a standard decoupling procedure based on 
the Hubbard-Stratonovich transformation in (|3.3| ), which introduce an auxiliary filed ip, one 
gets 



Z = C I V'tpVp exp 



1 /-i/r 



-J^ dr jdx{dM^ + [V-/^p 



+ TQip^ + Ipip^ - llP' 



1 

2uo 



(3.4) 



Using the fact that the field p is n-component we decompose the integral over ip into an 
n-dimensional gaussian integral, which can be performed easily, leading to 



Z = C I Vij) exp 



^^^-^TrHro + ^ + S^ + V^] 

4140 2 



(3.5) 



In the last expression we assumed that the field ip is time and space independent. For large 
n we use the saddle point method to evaluate the integral over ip in ( p.5| ). Finally, we obtain 
the free energy per particle T = — {T jY^ TilnZ in the momentum space 



J^ 



1 



4mo 



roY + 77 X! ^^ ^ ^^^^ 



V 



1 
2T 



+ Ifch 



(3.6a) 



and the saddle point equation 



T 



V 



k,m 



+ (27rmT)2 + |fc|'^' 



(3.6b) 



where, for convenience, we used the shifted parameter = ■?/; -|- tq instead of ip itself. 

Equations ( ^.6| ) are our starting point to explore the finite-size and quantum effects on 
the bulk critical behaviour of the model (p.l|) . Let us note that in the particular case of 
short-range forces a = 2, we recover the results of Ref. [l^ utilized to investigate the finite- 
size scaling when the quantum effects are absent. When the quantum effects are relevant, 
equation ( p.6b| ) was used in Ref. |T^ in order explore the quantum critical behaviour and 
to calculate the critical exponents. 

To extract the physics from model (|3.1|) , we calculate the susceptibility x- Iii the large 
n limit, we obtain 



X= d'^xi^ixMO)) = -. 



This is a particular case of the general correlator xik,uj) 
the dynamic susceptibility: 



X{k,UJn) = (0 



+ ujl + \kr 



(3.7) 
{(j){k,ujn)(p{-k,-ujn)), defining 



(3.8) 



From this expression one deduce that the anomalous critical exponent measuring the diver- 
gence of the correlation length at the quantum critical point is rj = 2 — a and the dynamic 
critical exponent is z = a/2. The denominator of the dynamic susceptibility ( p.8|) has poles 
in the complex fc-plane at A; = ± (— </> — (27mT)^) . The closest pole to the origin is the 
one corresponding to n = 0. This pole determines an exponential decay of the correlation 
functions. So the correlation length turns out to be 

e = <p-'^^- (3.9) 

91) one deduce the simple ratio between the critical exponents 7 and u to 



From (|3.7| ) and 
be equal to a. 

To investigate the bulk quantum critical behaviour of the model ( |3.1| ) at zero temperature 
we transform the sums in equations (|3.6| ) into integrals, i.e. 

d'^k 



tj: 



duo 

— and 

27r 



-E 



{2tiY 



to get 



and 



-Fn 



4Mn 



- r,Y + 



duo 

ro + uo I TT- 

Zn 



duo 



d'^k 

{27rY 



In 



+ cj2 + |fcr 



d'^k 



(3.10a) 



(3.10b) 



{27tY (f) + uo^ + Ikl"' 

These equations contains all the necessary information from which we can extract all what 
we need about the critical behaviour at zero temperature. Here we will pay attention to 
the scaling properties of the thermodynamic functions in the neighborhood of the quantum 
critical point given by 

uo f d'^k /2 



(2vr)' 



This integral is infrared convergent only for d > d^, where (i< = a/2 defines the lower 
critical dimension. It is easy to calculate all the critical exponents for the model we are 
considering here in the large n limit by expanding the right hand side of equation (|3.10b|) 
for small 0. In the expanded expression one can see a natural emergence of (i> = 3a/2 as 
the upper critical dimension. Here we will give some of these critical exponents: 

d — d^^ d — d^ 

In the remainder of this paper we will investigate the effects of finite temperature and/or 
spatial sizes on the bulk zero temperature critical behaviour. 



IV. FINITE TEMPERATURE EFFECTS ON THE QUANTUM CRITICAL 

BEHAVIOUR 

At finite temperature we obtain correction terms to the right hand side of equations 
( p.lOp . These are given by 



A^(T, 



for the free energy and 




(4.1a) 



AUT,, 




K, 1 m^ 



r 



OO -TV ^ _ j_ 

Err 2 
d_ 1 



(4.1b) 



for the saddle point equation. In expressions ( |4.1| ) we used the quantity k^^ = |(47r) 2r((i/2) 
and K^{x) is the MacDonald function (second modified Bessel function). 

Combining equations ( p.lOb| ) and ( [4.1b|) we get, in the bulk limit corresponding to the 
geometry oo'^ x L^, the following scaling form (for the saddle point equation) for dimensions 
d between the lower d^ and upper critical dimensions (i> 



, yi''-\(d' 



r 



1 




(4.2) 



Here the scaling variables are defined by yr = \/(p/T and Xr = T^~'^/^ (roc — tq). 

For fixed a, the solution, yo, of ( [4.2|) at the quantum critical point {xr = 0) is a universal 
number. An analytic solution of this equation cannot be obtained in the general case. Here 
we will consider some particular cases: 



'^{d-d<), t/-a/2<l. 



yo 



2 In 



2tt 



2 ' 

dy — d 
So- ' 



CT, 



(4.3) 



3(t/2 - rf < 1. 



At Xr = 0, equation ( [1.2|) can be solved numerically. The behaviour of the universal 
constant yo as a function of the reduced dimensionality d/cr is given in Fig. |I|. We see that 
yo depend upon the ratio d/a in a universal way for all values of a smaller than or equal to 
2. 



FIGURES 
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0.25 



reduced dimensionality d/a 

FIG. 1. Solution of equation ( |4.2| ) at x,- = as a function of the reduced space dimensionality 
of the system 

The phase diagram of the model (here we are speaking about the critical line and 
crossover lines) in the vicinity of the quantum critical point is determined by putting = 
in the saddle point equation ( ^■6b| ). In the vicinity of the quantum critical point ( |3.11[ ) the 
phase diagram is determined by 



TQc 



roc{T) = -k,T^''/'^-'^T 



a 



2d 



a 



1 c 



2d 



i. 



1 



(4.4) 



where ({x) is the Riemann zeta function. The crossover lines and the critical line are 
determined by It-qc — ''"ol ~ T^/'^^, The phase diagram is presented in Fig. ^. In the different 
regions of the phase diagram and for arbitrary dimensions the correlation length has different 
behaviors. 

For Xt — ^ C)0 in the interval cr/2 < d < a, the correlation length behaves like 



e 



T 



'^o - roc 



l/(d-a) 



(4.5) 



This expression shows that, when T -^ 0+, the correlation length goes to infinity as quan- 
tum effects becomes relevant. In this case there no quantum phase transition for finite 
temperature. For the second interval a < d < 3cr/2 one gets the behaviour 



e 



T 



l/id-a) 



ro-roc(T)^ 



(4.6) 



for ro very close but larger than the critical quantum parameter roc(T). Here we have a 
phase transition at finite temperature and the critical exponents are those of the model in 
the classical limit. They are: 



7 = 0" (d — (t) , u = (d — a) 
for To less than its critical value the correlation length is infinite. 



(4.7) 



9 



For Xr ^ —oo the correlation length is temperature independent. For arbitrary dimen- 
sion d^ < d < dy it is given by 

>i/{d<-d) 



The particular case d = 
simpler. Its solution leads to 



C ~ [roc - To) 
a is very simple to handle. In this case equation 



for the inverse correlation length. Here 

1 



/§ (xt) = 2 arcsinh 



exp 



a Xj 
'2E 



(4.8) 
2D get 

(4.9) 
(4.10) 



is a scaling function, which simplifies in some limiting cases: 



fd^) 



ax 

ZftfT 



1 l^ 1+V5 



exp 



I 2k„ ) 



X 



X 



X 



-00, 



(4.11) 



00. 



From equations (|4.11| ) one can transparently see the different behaviours of the correla- 
tion length ^(T) in three regions: 

a) renormalized classical region with exponential divergence as T — > 0. In this region 
the system display characteristics of the ordered ground state. The thermal fluctuations 
destroys long-range order at any finite temperature. 

b) quantum critical region with ^(T) ~ 2^-2/0- ^^^^ crossover lines T ~ |roc — tqI'^ • In 
this region the system "notices" that it is finite in the time direction. 

c) quantum disordered reg'^on with temperature independent correlation length as T ^ 0. 
In this region the system has a gap in the spectrum. 

The different regions are qualitatively shown in the phase diagram in Fig. ^. 

T 



Quantum 
critical 



Renormalized 
classical 



Quantum 
disordered 



FIG. 2. Qualitative phase diagram of the model under consideration in the vicinity of the 
quantum critical point for d = a. The dashed lines are crossover lines. The tick line shows the 
region, where the system is ordered at T = 0. 
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V. QUANTUM CRITICAL FINITE-SIZE SCALING 



As the temperature is set to zero, equation (|3.6| ) turns into 



1 



^L = :^Y.^J^+\kV-^^{^- rof (5.1a) 



for the free energy and 






for the saddle point equation. 

Exploring the finite-size scaling in the quantum limit (T = 0) turns out to be a difficult 
task because of the presence of the term Ifcl*^ in the spectrum of the model. In other words 



the problem is how to handle the terms \/(p+ \k\^ in the expressions for the free energy 
and in that of the saddle point equation ( p.l| ). The solution to this problem was given in 
Ref. [^, where two important identities facilitating the analysis of the finite-size scaling 
were presented (see also Appendix ^. 

For the (i-dimensional system with spatial geometry L'^^'^ x oo'^ at zero temperature and 
periodic boundary conditions imposed along the {d — d!) finite size directions with linear size 
L of the system, the corrections to equations (|5.1|) are given by 



and 



2 J (2'Kf U^ Ufga 



where for convenience we introduced the following functions 
and 

Here the primed summation over the vector Z is (d — (i')-dimensional and the prime indi- 
cates that the term corresponding to Z = is excluded. In the last equation we used the 



function lEO 



1 ^ r(A; + l/2) x^ 
QaA^) = ^=22 r/^,. , m TT (5.2e) 



;^o r {ak + (3) k\ 



Some properties of the functions Qa,i3 and A|(y) are discussed in Appendices and 0, 
respectively. 
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A. The finite-size shifted critical quantum parameter 

It is well known from finite-size scaling theory that the critical value of the parameter 
driving the phase transition is shifted due the effects of the finite sizes in the system. The 
aim of this subsection is to evaluate the distance over which the critical quantum parameter 
roc of the bulk system is shifted. For our concrete model this is obtained by substituting the 
parameter by zero. The result is 

_ Up Tid/2 - a/A) ^, (\l\L\-'^''/' 



and the {d — d') dimensional sum in the right hand side of equation ( |5.3|) is convergent for 
d' > ci<. 

In the opposite case (i< > d', the r.h.s of ( |5.3|) is divergent. Nevertheless the sum in 
equation ( p. 3D can be expressed in terms of the Epstein zeta function, which is a gener- 
alization of the Rieman zeta function. The resulting Epstein function can be analytically 
continued beyond its domain of convergence to give a physical meaning to equation ( ^.3|) as 
well. In this case the shifted "pseudocritical" quantum parameter roc{L) corresponds to the 
center of the rounding of the singularities of the thermodynamic functions, taking place in 
the thermodynamic limit. This point has been investigated in details in the framework of 
the finite-size shift of the critical temperature for the spherical model in Ref. pl |. 



Let us also mention that the distance over which the critical quantum parameter is 
shifted can be also expressed as 

roc - ro{L) = 2(2vr)<^/2 r(a/4)^^'^'-' ^^'^^ 

where Cd,d',a is a the Madelung type constant (cf. Appendix 0). 

One can see that the shifted critical quantum parameter rQc{L) is lower than its bulk 
bulk critical value tqc for the different values of d, d' and a (which is the "normal case" , see 
Ref. 10), while the pseudocritical rQc{L) is upper than the bulk critical quantum parameter. 
However for the boundary case when d' — >■ (i< we find that the shift is infinite. This may 
be explained with the aid of the behaviour of the Epstein zeta function at its pole d' = d^. 
The shift in this case is 6rQ ~ {d' — d^)~^ and the appearance of =Foo is clear. 



B. Finite-size scaling at zero temperature 

In the neighborhood of the quantum critical point, it is possible to write equations ( |5.2| ) in 
the scaling forms (for dimensions between the lower d^ and the upper dy critical dimensions) 



^i - ^0 = L-''-' 



-DlA^"" + l^^^LVL- ^"^riVL) 



(5.5a) 



for the singular part of the free energy and 
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Mo 



XL = ^Dd^„yl ^+MoA|^(yz.) 



(5.5b) 



for the saddle point equation, whose solution is related to the correlation length. Here the 
scaling variables are given by xl = {tqc — ro)L''~'^< and yi = L^cj), and 



D, 



d.o 



-M^i 



1.^ 






(5.5c) 



From equations ( |5.5| ) one can see that the singular part of the free energy and the 
correlation length are universal scaling functions of the variable xl i.e. 



J-',^ = J-" — J-'o 



-d—z 



fA^L) 



and 



e = Lf^ixL) 



(5.6a) 



(5.6b) 



At the quantum critical point i.e. xl = the critical amplitudes fj^{0) and /^(O) are 
depending upon the geometry of the system and the range of the interaction. 

i) In the region corresponding to X/, — > cxd (in other words if the quantum parameter rg 
is less but very close to its critical value tqc), we obtain 



e 



Joc-ro \''<' 



^Hr 



UqD, 



d',(7 



(5.7) 



This shows how the correlation length tends to infinity as the size of the system becomes 
larger. In the bulk system we recover the fact that the correlation length is infinite in the 
ordered phase in the large n limit. 

ii) In the case when the quantum parameter tq coincides with its critical value roc, the 
value of the scaling function /^(O) determines the critical amplitudes at the critical point. 
The value of these critical amphtudes depends upon the dimension d of the system, the 
number of infinite sizes d! and the range of the interaction a. It can be evaluated analytically 
in the vicinity of the borders (determined by the critical dimensions) where the scaling is 
valid. The scaling variable vanishes and the function Ag(?/i), in equation ( |5.5b| ), can be 
replaced by its asymptotic form for small argument ( P4a|) . As a solution for the obtained 
equation we find 



L 



(T/2-d' 



.(4^)^(d-f)r(f)D^,, , 



I V(4^)T(3f_d)r(3^)D, 



3(T/2-d' 



rf-f <1, 



f -t/<l. 



(5.8) 



In some particular cases when the function Qa,ii{x) gets simple, namely in the cases a = 1 
and cr = 2 (see Appendix 0), and for some special cases of the dimensions d and d' it is 
possible to solve equation ( |5.5b| ) numerically. Then one gets 
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i 

L 



( 0.624798 for d = 1, d' = 0, a = 1 

1.511955 for d = 2, d' = 0, a = 2 

^ 2 In (^) for d = 2, d' = 1, a = 2 



(5.9) 



in) The last case we consider here is the one corresponding to the values of the quantum 
critical parameter smaller that the critical value i.e. xl -^ — cxd. The correlation length is L 
independent and it is given by 



UoDd,a 

2{roc-ro)^ 



(5.10) 



VI. INTERPLAY BETWEEN QUANTUM AND FINITE-SIZE EFFECTS 



In this section we consider the case of a finite system at low temperatures. In other words 
we will investigate a system confined in the finite geometry of the general form L'^"'^ x cxd*^ x 
L^ (here we will limit our discussion only to the case d' < d^). In this geometry, in addition 
to the correction terms to equation ( ^.10b| ) given in ( [4.1b| ) and ( |5.2dD we have to add also a 
correction term which accounts for the combined effects of finite sizes and finite temperature. 
This is given by 



T (0, L, T) = j^ J2J2' -T exp(-2<; 



2^2 



where 



(27rj ^ m I 

x|Zr''$rf/2-i,.(;2L-'^ir" 



$,,.(y)= / rfxx'^+V,(x)e-^^"' 

JO 






(6.1a) 



(6.1b) 



was introduced in Ref. [^. In equation (|6.1b| ), J,y{x) stands for the Bessel function. 

The general equation obtained by putting all together equations (|3.10b| ), ([4.1b|) , (|5.2d|) 
and ( |6.1| ) can be written in a scaling form whose solution gives the correlation length as a 
function depending upon two scaling variables. This has the general form 



e = Lf^ {xL, LT^'^) = T-'/^U (a:., LT^ 



1-^ 



(6.2) 



Actually we see that there will be a competition between the finite sizes and quantum 
effects depending upon the quantity LT^I^ as we will see latter. The scaling form (|6.2|) is in 
agreement with predictions of Section 0. 

The solution of the saddle point equation for a system confined to a general geometry 
is very difficult to obtain in an explicit form. Even in the two limiting cases of "low- 
temperature" {LT^I^ ^ 1) and "very low-temperature" {LT^I^ -C 1) the asymptotics of 
the general equation are very complicated. Nevertheless these limits provide some useful 
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information about the behaviour of the system in the first or the latter case. Let us mention 
that the ensuing mathematical equations simplifies drastically in the case of short-range 
forces i.e. a = 2. In this particular case the analysis of the saddle point equation is identical 
to the one presented in the framework of the quantum rotor spherical model (for details see 
Ref. [jl^)- There is a case when the saddle point equation takes more simple form. This is 
namely the particular case of long-range interaction corresponding to cr = 1. Even in this 
case the analysis of the critical behavior is very complicated. 

Before doing this let us present the expressions for the asymptotic forms corresponding 
to the limiting cases when the finite temperature effects dominate the finite-size scaling and 
vice versa. 

A. Low-temperature regime LT^'''^ ^ 1 

In this regime the finite temperature corrections lead those coming from the finite size 
of the system. The scaling form of the saddle point equation is given by 

where we used the functions 




/C2(z/; y) = 2 £ imy)-''K,i2my) (6.3b) 



m=l 

and 



r- f \ ^' f°° , x^^JA2ylx) 

M-;y) = E / dx ^ ^^ . (6.3c) 



l + x"" (yl) 



The last function has been proposed in Ref. [^, where the finite-size scaling properties of 
the spherical model were discussed. The analytical properties of this function are considered 
in the same reference. 

The function /C2(i^; y) introduced by equations ( |6.3D is exponentially decreasing for large 
values of its argument y (for fixed finite u). The second function, i.e fC„{u,y), is decaying 
exponentially only in the case of short-range interaction a = 2. For the case of long-range 
interaction corresponding to a < 2 its asymptotic form decreases by a power-law. The two 
functions are identical in the case a = 2. 

Let us recall that the temperature enters in the right hand side of equation ( |6.3a{ ) trough 



the relation y^ = y/(f)/T. The last term appearing in the scaling form (|6.3a| ) is a correction 



to the bulk system at low temperature resulting from finite-size effects. The temperature 
independent part of this term is nothing but the finite-size corrections to the bulk system in 
its classical limit i.e. the case when the quantum effects are irrelevant. Now, in the case of 
short-range interaction {a = 2) this equation has been analyzed previously [1^. In this case 



we get exponential corrections. In the case of long-range interactions we have a power-law 

15 



corrections. This seems to be a general characteristic for systems with long range interaction 
(for the case of the spherical model see references |23,2^). 



In the particular case of the two dimensional system with short-range interaction d = 
a = 2, we get the solution 

r' ^9T + {2- d')^^^j exp i-TLO) , (6.4) 

where 6 = .962424 is a universal constant. This universal form is agrees with the predictions 
announced in Section 0. 

Let us quote here the result for the shift of critical quantum parameter from its bulk 
value. It is 

Toe - ro.mT) = ^uok,T(''^/-'^T (^ " l) C (^ - l) 



2T{A7TY/^r{a/2)^i 

The first term in this equation is the finite temperature shift of the critical quantum param- 
eter and the second term is nothing but the finite-size shift in the classical limit (i.e. when 
the quantum fiuctuations are negligible in the system) divided by the temperature. 

B. Very low-temperature regime LT^'^ ^ 1 

Here the finite temperature correction are negligible in comparison with those coming 
from the finite-size effects. In this case the scaling form of the saddle point equation turns 
into 

XL = Y^d,ayL + Uo\ [yL) 

From this equation we can see that the finite temperature corrections to the finite system 
at zero temperature are coming mainly from the d' infinite dimensions. In this case the 
temperature corrections are exponential. In the particular case of two dimensional system 
with short-range interaction d = a = 2 confined to a strip geometry {d' = 1) at the quantum 
critical point, we get 

^ , e [2^ [l ( 9 \ , „, 

In the case when the system is confined to a block geometry {d' = 0), we have 



,1^1 
£"- — + - 
^ L L 



2fi 2 Y ^ 
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-1 



exp ( --^ ) , (6.8) 



where fi = 1.511955 is a universal constant. Both results ( |6.7|) and ( |6.8|) are in agreement 
with the scaling predictions of Section |T|. 

The shift of the critical quantum parameter in this case is given by 

^0 T{dl2 - a I A) , (WA-'^""" 

Toe roel^,JJ- 2 (4^),/2r(^/4)Z. 1^ 2 ) 



d'ri.{2d'/a~l)-p /^ _ ^\ ^ /^ _ A ^ 



--kd'L''~'''T^"''l"~'^V (=^-l)C(=^-l). (6.9) 



The first term of the right hand side of equation ( |6.9D is discussed in subsection |V A 



The second term is the correction to the finite-size shift coming from the temperature. It 
is just the finite temperature shift of a d! dimensional system multiplied by the volume of a 
[d — (i')-finite hypercube with linear size L. 

C. Short-range case a = 2 

The equation for the saddle point equation reads 



Tqc - To 



""O T^ I ^ ~ d\ ^('rf_i)/2 



(47r)^ V 2 



(47r) '2' ^^j J01/2 (^2/^2 + ^2^; 



This is the simplest form one can get for the saddle point equation. In this case the dynamic 
critical exponent z, which measures the anisotropy is 2; = 1. The system is symmetric with 
respect of the change L ^^ T^^. 

The general analysis of equation ( |6.10[ ) follows the one presented in the framework of the 



quantum rotors model. Here we will not discuss this point since one can see the details in 
Ref. [ITTi. 



VII. SUMMARY 

The 0{n) vector (p^ model is extensively used in the analysis of the critical phenomena 
because of its direct relevance to the physical reality. In the limit ri — > 00 it adds the property 
of exact solvabilty at any dimension. This makes it very attractive for the exploration of the 
the scaling properties of quantum critical phenomena as well as finite-size scaling theory. 

In the paper we present investigations on the finite-size scaling of the yj^-model in the 
vicinity of its quantum critical point. The striking characteristic of the model is the presence 
of long-range interaction decaying at large distances r with a power law as r~'^~'^ . We 
considered the model confined to the general geometry of the form L'^"'^ x 00 '^ x L^, where 
L is the spatial size of the system, L^ ~ T~^/^ and z is the dynamic critical exponent. 

The results are obtained by considering the temperature which governs the crossover be- 
tween the classical and the quantum critical behaviors as an additional temporal dimension. 
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A detailed investigation of the alteration of the zero temperature critical behaviour of the 
model due to the finite temperature is performed. For dimensions (i< < d < (i>, we studied 
the critical behaviour of the bulk model in the tree different regions: classical renormalized, 
quantum critical and quantum disordered with different behaviors of the correlation length 
as function of the temperature. The behavior of the correlation length in these different 
regions as a function of the dimensionality is calculated. Also, some critical amplitudes are 
evaluated. 

The large L behavior and the scaling forms accounting for the finite size effects is derived 
for the free energy and the saddle point equation. For the finite-size shift of the the critical 
quantum parameter, we found tq^ — '"o ~ L^^^'^ in accordance with the general postulates 
of finite-size scaling. For some particular cases the behavior of the correlation length and 
some critical amplitudes are evaluated. 

The study of the general case when the system is confined to the general geometry 
L'^^^ X oo'^ X L^ i.e. when the temperature as well as the sizes of the system are finite 
turns out to be a very difficult task because of the high anisotropy of the system due to the 
parameter a. Nevertheless one can make interesting deductions in some limiting cases: (i) in 
the low-temperature regime [LT'^^^ ^ 1), the finite-size corrections to the low-temperature 
behavior are exponentially small in the case of short-range interaction and are decreasing 
with a power law in the case of long-range interaction, (ii) in the very low-temperature 
regime {LT^^^ <^ 1), however, the finite temperature correction to the finite-size behaviour 
are always exponentially small. 

Here, we confined our investigations on the static critical properties of the model, be- 
lieving that the spherical limit (n = oo) provides a useful tool for studying quantum critical 
phenomena in dimensions d > a/2. The dynamic properties are not well described in this 
limit and require loop corrections (see e.g. reference ^5\ in the case of the nonlinear quantum 
sigma model). 

Though derived for the special case of the ip"^ model with long-range interaction in the 
large n limit, the obtained here results are expected to hold also for many other cases. For 
example, recently, a model suitable to handle the joint description of classical and quantum 
fiuctuations in an exact manner, was considered in a number of publications [p0|,p6|-|30 



This model is a modification of the y9"^-lattice model used extensively in the investigation 
of the critical behaviour of the the structural phase transitions [^, in the spirit of the 
self-consistent phonon approximation method |2^ . 
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APPENDIX A: SOME PROPERTIES OF THE FUNCTION GaA^) 



The functions Ga,i3{z) were introduced in Ref. |2^ in order to investigate the zero tem- 
perature finite-size scaling of an anharmonic crystal with long-range interaction at zero 
temperature. They are entire function series of finite order of growth defined by 



^'^At) = ^t VTTT^VV « > 0, /3 > 0. (Al) 

V71" ^ r {ak + (3) k\ 

One of the most striking properties of this function is that it obeys the following identity 

-^ = / dxe-'^x^-^g^^p (-x°^) , (A2) 

\/l + z Jo 

which is obtained by means of term- by-term integration of the series ( |AT|) . The identity (|A2|) 
lies in the basis of the mathematical investigation of finite-size scaling in quantum systems 
with long-range interaction. 

In some particular cases the functions Qa^piz) reduce to known functions. Here we will 
give some examples which are of interest for us in this paper: 

Gi,i/2{z) = ^, (A3a) 



Gi/2,i/4i-z) = -^U (^, i, z^) , for ;. > (A3b) 



4^/^ V4'2 

and 

1 ../I 1 



6^iA3/4(-^) = -^U (^-, -,z'], for z > 0, (A3c) 



where U{a, b, z) is the confluent Hypergeometric function |^| 



If we set in the identity (|A2| ) z = y ", y > 0, and x = ty, we will obtain the Laplace 
transform 



'l + yay/2 =l dxe-yH^-'g^A-n, (A4) 



,«/2-/3 ^oo 

yC 

from which we derive a new identity by setting (3 = a/2: 

(l + ^-)-V2=/ dxe-''x''/^-'g,„/2{-xn- (A5) 

Jo 

With the aid of the last equation we obtained the large L asymptotic behaviour (|5.2b[) from 



equation (|5.1b|) . 

The integral representation of the functions ga,i3{z) is obtained with the aid of the Han- 
kel's integral for the inverse of the gamma function 

^ ^ ^ e^'u-'dz, (A6) 



T{z) 2Txi Jc 

where the integration contour C is a loop which starts and ends at x = — oo and encircles 
the origin in the positive sense: — vr < argz < vr on C [^. This enables to get the result 



Gasiz) = 7-- / dv—. m^- (A7) 

'^^ ' 2-KlJc (W°-Z)V2 ^ ' 

19 



The last result is valid only for a < 1. The integrand in ( |A7|) has a branch-point at 
f = 0. A more complete and detailed analysis of the function Qa,i3{z) will be presented 
in a subsequent paper. This integral representation is helpful for obtaining the asymptotic 
behaviour for z ^ oo. 

In the following we will investigate the asymptotic behaviour of the functions Ga,f3{—z) 
for real argument z. This may performed by the use of the series 

1 _ 1 ^ r(fc + i/2) / ty ^ 1 g r(A: + i/2) / ty ^^g^ 



Vx + t \f2x^^ k\ \ xj V2x ^p k\ \ X 

ioY X ^ t . The second part of the right hand side of the last equation is nothing but 
the expansion in series of the hypergeometric function 2Fi{a,b;c,x) multiplied by some 
coefficients. So by simple rearrangement equation ( |A8|) takes the form 

1 1 '^r(A; + l/2) f t^"" 



y/x + t \/2x rr'n k\ 



X 



1 / ;z\P+i r(n + 3/2) / 3 z\ ,, , 



Using the integral representation (^) and the identity ( [A9| ) we get the asymptotics (for 
p>l) 

^"'"^ ' to k\^ TiP-aik + 1/2)) 

+0 (x-P-^/^) , X -^ +00. (AlO) 



In the particular case (3 = a/2, the last equation reduces to 

-3/2 



&./.(-) --5^ (All) 



APPENDIX B: ASYMPTOTICS OF THE FUNCTION As:{y) 

Here we present the asymptotic behaviours of the functions A^{y), given in ( ^.2d|) for 
small and large y. These functions are defined by 

xx''"-"-'g,/,,,/t {-x"-'y) . (Bl) 

With the aid of the Jacobi identity for a d-dimensional lattice sum 

Ee-'^=g)'^'Ee-^'''/^ (B2) 
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we transform expression ( pT|) into 



1 TT 2 



2 27r 2 7o 2'4 



X2 



{2TTy 



j:'' 



^xP 



xj 



(B3) 



In order to be able to get a reasonable expression for the integral in the last equation we 
have to avoid the divergence in the square brackets. To this end we add and subtract from 
the function Qa,i3{x) its small asymptotic behaviour to the first order (see Appendix 0), 
which enables us to write (after some algebra) 



A^iy) = ]^Da',y'"~"' - ^D.^.y"/''-^/^ + W,,,>Ay) 



+ 



vr' 



d'/2 



(27r)'^ r(CT/4) 



Cd,d',a, 



(B4a) 



where we used the notations 



>^..',-&r''-""-^''-'. 



-xi-' 



(2^) I 

Qa/2,(7/4 



X 



x'^/^y' 



(27r)-y r(a/4) 



(B4b) 



Cdd'a = E' r dxx''"-''"-\-^'" - -n^'-"'^'^ r dxx""-"''^-' 

' ' I Jo JQ 

= lim|y'rf---,A/2^|Z|^V2-./4 

'-Hi V4 2' y" 






(B4c) 



Here T{a,x) is the incomplete gamma function. The expression in equation ( [B4cD is a 
generalization of the Madelung type constant. It is y independent. Indeed this equation 
defines the finite-size shift of the critical quantum parameter tqc at zero temperature. It is 
easy (following Ref. [|1^) to show that equations ( P4c| ) and ( |5.3| ) are equivalent. 

For small y the asymptotic behaviour of A^(y) is given by the first term in the right 
hand side of equation ( |B4a| ) . 

For large y the asymptotics the function A^(?/) are obtained by substituting the large x 
behaviour of the functions Qa,i3{x) from equation (|A1U ) in the definition ([Bl|). After some 
calculations one ends up with 



My) 



-3/2 



T{-a/2) \ 2 )^\7r\l 



d+cr 
2 



(B5) 
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